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ABSTRACT 

We perform phase-resolved spectroscopy of the accreting millisecond pulsar, SAX 
J1808. 4-3658, during the slow-decay phase of the 2002 outburst. Simple phenomeno- 
logical fits to RXTE PCA data reveal a pulsation in the iron line at the spin frequency 
of the neutron star. However, fitting more complex spectral models reveals a degen- 
eracy between iron-line pulsations and changes in the underlying hotspot blackbody 
temperature with phase. By comparing with the variations in reflection continuum, 
which are much weaker than the iron line variations, we infer that the iron-line is not 
pulsed. The observed spectral variations can be explained by variations in blackbody 
temperature associated with rotational Doppler shifts at the neutron star surface. By 
allowing blackbody temperature to vary in this way, we also find a larger phase-shift 
between the pulsations in the Comptoniscd and blackbody components than has been 
seen in previous work. The phase-shift between the pulsation in the blackbody tem- 
perature and normalisation is consistent with a simple model where the Doppler shift 
is maximised at the limb of the neutron star, ~ 90° prior to maximisation of the 
hot-spot projected area. 
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1 INTRODUCTION 

SAX J1808. 4-3658 was the first accreting millisecond X-ray- 
pulsar (AMXP) to be discovered (jWijnands & van der Klis 
19981. Lying some 3.5 kpc distant (Galloway & Cumming 



2006), the neutron star spins with a frequency of ~401 Hz 



and is thought to be orbiting a low-mass brown-dwarf com- 
panion (jBildsten fc Chakrabarty|2001[ ) with a 2.01 hr orbital 
period ( |Chakrabarty fc Morgan|1998[ ). The established pic- 
ture is that such AMXPs are recycled to millisecond periods 
from more slowly rotating neutron stars, due to the trans- 
fer of angular momentum from the accretion flow provided 
by the companion (Bhattach arya fc van den HeuveT||1991[ ). 
In this sense, the AMXPs bridge the gap between low mass 
X-ray binaries (LMXB) and rotation powered millisecond 
pulsars. 

The accretion flow in AMXPs differs from that in black 
hole X-ray binaries (BHXRB) due to the stellar magnetic 
field and the solid neutron star s urface. The magnetic field 
of SAX J1808. 4-3658 is ~ 10 s G ( |Hartman et al.|2008| > and 
this is sufficient to channel accreting material approaching 
within a few neutron star radii along field lines towards the 
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magnetic poles. The funnelled gas passes through an accre- 
tion shock before falling on to the magnetic poles. X-ray ob- 
servations of AMXPs made with the Rossi X-ray Timing Ex- 
plorer (RXTE) and XMM-Newton (XMM) have discovered 
distinct spectral components originating from these sources. 
The accretion disc and neutron star surface contribute to 
a soft component up to ~10 keV, and Comptonisation in 
the shock region (most likely of seed photons from the im- 
pact hot-spot) produces a harder component extending to 
~100 keV ( jGierlinski, Done, fc Barret|2002[ ). Recent XMM- 
Newton observations of SAX J1808. 4-3658 have shown two 
distinct blackbodies at soft energies plus a hard component, 
interpreted as an accretion disc, a neutron star hot-spot and 
an accretion shock ( Patruno et al.|2 009 Pap itto et al.|2009 1. 

In an analogous manner to reflection in BHXRB 



( |Fabian et al. 1989), it is thought that Comptonised pho- 
tons reflect off the thin optically thick accretion disc pro- 
ducing a Compton reflection hump and fluorescent iron Ka 
line emission at ~6.4 keV (see Cackett et al.|2009a for a re- 
view of broadened iron lines in neutron star LMXBs). The 
shape of the line provides information about the distance 
of the line-forming region, since reflection occurring close 
to the compact object produces a characteristic asymmet- 
ric profile due to Doppler-broadening and gravitational red- 
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shift. Observations of SAX J1808. 4-3658 have revealed evi- 



dence suggestive of such a broadened iron line (Cackett et 
al.|2009b||Papitto et aL|2 009). Iron Ka lines are also of par- 
ticular importance because they place a constraint on the 
magnetospheric radius. It is, as yet, unclear how the mag- 
netic field of a magnetised star truncates the disc, or what 
effect the star's rotation has on the disc inner edge. 

The X-ray emission from AMXPs is pulsed due to mis- 
alignment of the spin and magnetic poles. As the hot-spot 
and magnetically channelled accretion column rotate around 
the spin axis, the projected area and the magnitude of the 
Doppler shift changes, causing emission to be pulsed at the 
spin frequency. It is therefore not unreasonable in this sce- 
nario to imagine that changes in the line shape and/or en- 
ergy with phase might also occur. For instance, if the accre- 
tion column illuminates significantly different parts of the 
disc with phase, reflection occurring from approaching and 
receding parts of the disc should display a change in line 
energy with phase. Exploring the variation of the line prop- 
erties with phase, as well as the phase dependent correlation 
between the line equivalent width (EW) (or normalisation) 
and the reflection hump, might be a valuable window on the 
formation of this spectral feature. The iron Ka feature in 
SAX J1808. 4-3658 has been known about for many years 
(Gicrlinski, Done, & Barret 2002), and we now know that 
this feature might be relativistically broadened (Cackett et 
al.||2009b| |Papitto et al.|| 2009[ ). Therefore, it is of interest 



whether this feature is also pulsed at some level, which might 
be expected if the illuminating continuum is also pulsed as 
seen from the inner disc (which might be expected in certain 
geometries, e.g. where a part of the disc is shielded by the 
neutron star itself). 

Using the correct orbital ephemeris and referencing the 
photon arrival times to the Solar System barycenter, it is 
possible to define an arbitrary phase-zero and 'fold' all of 
the incoming photons in to phase-bins representative of a 
full rotation of the neutron star. By examining the energy 
spectra of photons from individual phase bins, one can per- 
form phase-resolved-spectroscopy. A phase-resolved analy- 
sis was performed of the 1998 outburst ( |Gierlinski, Done^] 
& Barret 2002 ) , and the self-consistent continuum reflection 
and iron-line were found to be consistent with having a con- 
stant value as a function of phase. In this work, we produce 
phase-resolved spectra of the best-sampled and best-studied 
2002 outburst of SAX J1808. 4-3658 as observed by RXTE 
over the 'slow-decay' of the outburst ( |Hartman et al.||2008[ ) 
where pulse profiles are most stable and we can investigate 
the phase-resolved spectrum with high signal-to-noise. We 
shall explore changes in the iron-line properties with phase, 
and investigate the viability of the currently accepted pic- 
ture of the broadened iron-line emission in AMXPs. 



2 OBSERVATIONS AND DATA REDUCTION 

We have concentrated on RXTE observations of SAX 
J1808.4-3658 from the slow-decay phase (MJD 52564-52574) 
of the 2002 outburst where the pulse profile is the most sta- 
ble and the data is best sampled (in date order OBS ID 
70080-01-01-02 to 70 080-01-03-03 inclusive ). The reader is 
referred to Figure 3 of |Hartman et a l. (2008) for a light curve 
of the outburst illustrating the data segment we used. For 



simplicity, only data taken with the E_125US_64M_0_lS 
event mode of the RXTE Proportional Counter Array 
(PCA) were used in this work. Although during these obser- 
vations different PCA detectors were switched on, we have 
only analysed data from PCU 2 which was always on, in or- 
der to ensure consistency. Gain and offset corrections, which 
attempt to match the energy scales in different PCU detec- 
tors, assign zero response in different channels of different 



PCUs (Jahoda et al. 19961. Combining different detectors 
involves weighting different detector responses, and this in- 
troduces unmodelled scatter in to the data. In total, 64 event 
files were analysed with a total exposure of 167 ksec, which 
provided sufficient signal to noise for phase-resolved spec- 
troscopy using 16 phase bins. 

The photon arrival times in each event file were 
corrected to the Solar system barycenter with the ftool 
faxbary. The photon phases were determined by using a 
high-order polynomial to "whiten" the phase residuals and 
increase the signal to noise of the pulsations. Additionally, 
the arrival time of photons was corrected for the orbital mo- 
tion of the binary using a first-order approximation as out- 
lined, for example, in |Papitto et al.| ( |2005[ ). A folded pulse 
profile was created for each separate event file, and fitted 
with a simple constant offset plus a sinusoid. During the 
slow-decay, pulse profiles have been shown to be stable and 
nearly sinusoidal ( Ibragimov & Poutanen 2009 Hartman et 
al]|2008 1. By fitting each event file with the best fitting si- 
nusoidal function, an offset could be determined for each 
event file to accurately align the different files in phase for 
subsequent summation. Custom code was written to select 
the photons from PCU2 only and produce separate FITS 
files containing events from each of the 16 different phase 
bins for each event file. The ftool SEEXTRCT was then used 
to produce PHA spectra for each phase bin. 

Background files were created using the ftool PCABACK- 
est, and extracted using SAEXTRCT over contemporaneous 
good time intervals created using the ftool maketime. The 
ftool mathpha was used to sum pha files from each respec- 
tive phase bin and sum the background spectra. A weighted 
PCU2 response file was created (due to any variations in 
the response over time) from all observations using PCARSP 
and ADDRMF. Care was taken to adjust the exposure of the 
phase-resolved spectra to reflect the fact that each spec- 
trum only contained l/16th of the total exposure. Phase 
averaged PCU2 spectra were also obtained from the same 
observations. 

The pre-processed High Energy X-ray Timing Experi- 
ment (HEXTE) cluster A spectra from the same observa- 
tions were summed together to provide some constraints on 
the phase-averaged power-law spectral component at higher 
energies. 

Throughout this work, spectral fitting was performed 
using XSPEC version 12.5.1 ( |Arnaud|p9lil ) with a 0.5% 
systematic error for PCA and HEXTE data during fitting. 
It should be noted that the PCU 2 data has zero response 
in channel 10, but these counts are redistributed to other 
channels, hence the gap in the data in Figure [2] 
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3 ANALYSIS AND RESULTS 
3.1 Phase-averaged spectrum 

Phase-averaged spectra of SAX J1808. 4-3658 from the slow- 
decay phase of the 2002 outburst have been fitted before 
| |Ibragim ov fc Poutanen 20091. It is therefore known that 



a simple power-law fit shows residuals that are consistent 
with Compton reflection and blackbody hotspot emission. In 
both the 1998 and 2002 outbursts, thermal Comptonisation 
models have typically shown photon indices of ~ 1.8 with 
blackbody temperatures a round 0.7 keV and hot-spot areas 
in the range 20-110 km 2 (|Gierlihski, Done, fc Barret||2002[ 



Ibragimov fc Poutan en 2009; Pout anen fc Gierlihskif2 003 ) 

As well as fitting PCA data over 3.0 to 25.0 keV, the 
HEXTE spectra are also fitted simultaneously over the range 
15.0-150 keV to constrain the power-law continuum and cut- 
off energy. The continuum was best-m odelled using the ther- 
mal Comptonisation model nthcomp (Zycki, Done, & Smith 



19991, and reflection was modelled using pexriv ( Magdziarz 



& Zdziarski 1995 1 and GAUSSIAN components convolved with 



rdblur. Only the reflection component of the pexriv model 
was included. A thermal Comptonisation model is chosen 
since the hot electrons in the accretion shock are thought 
to Compton up-scatter the seed photons from the hot-spot, 
whilst the blurred Gaussian is chosen to model the iron- 
line feature. All parameters were tied between the two in- 
struments apart from a normalising constant to allow for a 
calibration offset. The hotspot emission is modelled using 
a BBODYRAD component. The bbodyrad model is a sin- 
gle temperature blackbody spectrum parameterised by the 
temperature and the normalisation Rk m /D 2 , where Rkm 
is the source radius in km and Dio is the distance to the 
source in units of 10 kpc. The nthcomp model is parame- 
terised by the asymptotic power-law photon index, the elec- 
tron temperature, the seed photon temperature and the nor- 
malisation (the photon flux defined at 1 keV in units of pho- 
tons cm _2 s _1 keV _1 ). The pexriv model describes an expo- 
nentially cut-off power-law spectrum reflected from ionised 
material and is parameterised by the power-law photon in- 
dex, the power-law cut-off energy, the disc temperature, the 
disc ionization parameter ( Done et al.|1992~ l, the normalisa- 
tion (photon flux at 1 keV) of the power-law only (without 
reflection) and a dimensionless reflection scaling factor, rd- 
blur is a convolution model taking account of the blurring 
due to relativistic effects from an accretion disc around a 
non-rotating black hole, and is parameterised by the power- 
law index of emissivity, the disc inner radius, the disc outer 
radius and the inclination. Full parameters of the fit are 
shown in Table [l] and the unfolded spectrum is included in 
Figure [T] The spectral fits reflect the average parameters 
over the 167 ksec of data we analysed from the slow-decay 
phase of the outburst. The parameters of this 'average' fit 
are consistent with previous work, though it should be noted 
that as the mass accretion rate drops throughout the out- 
burst, the disc inner edge likely recedes (Ibragim ov fc Pouta-| 



20091, causing changes in the fit parameters. The fits 



were insensitive to £, the ionisation parameter, and the disc 
temperature so these were frozen at values of 100 erg cm 
s _1 and 10 6 K respectively. The low value of £ is consistent 
with the line energy which was constrained to lie between 6.4 
and 6.9 keV and pegged at the lower limit during fitting. In 
order to avoid the pexriv power-law cut-off energy becom- 



Table 1. Best fitting parameters for the phase-averaged 
consxphabsx (nthcomp+bbodyrad+rdblurx (pexriv+gaussian) 
model used in section [3JJ 



Parameter 



Value 



BB kT (keV) 
BB norm 

r 

nthComp KT C (keV) 

nthComp norm 

foldE(keV) 

PEXRIVnorm 
PEXPJV rcfl 

PEXRIV 5 
HEXTE C onst 

Rin (RG) 

Emissivity 
Fe line E (keV) 
Fe line width (keV) 
Fe norm (10 -4 ) 
X 2 / d.o.f 



0.676 
130.89 
1.96 



+0.02 
-0.01 
10.12 
12.87 
+0.03 
-0.03 
+ 14.08 



45.53 

051 + - 004 
u.uoi_ Q04 

136 5Q+ 42 24 

0.034 (Frozen) 
o oq+0.93 

100 (Frozen) 
ISAIAH 



-3.30 
6.40 



+1.1 

.7* 
+0.08 



-0.00* 
O.Ol(Frozen) 

6 H+ 2 - 70 

97.18 / 62 



Inclination was fixed at 60°, Rout fixed at 1000 Rq, nH fixed at 
0.113 X 10 22 and solar abundances assumed. Upper/lower error 
limits labelled with an asterisk show that the fit parameters 
were pegged, at times indicating that the fit was not 
particularly sensitive to this parameter. 



ing unphysical, this parameter was constrained to be equal 
to three times the value of the electron temperature. Free- 
ing the seed photon temperature for Comptonisation from 
the blackbody temperature had no statistically significant 
effect, consistent with the notion that the seed photons for 
Comptonisation come from the hotspot emission. It should 
also be noted that the pexriv normalisation and reflection 
fraction are degenerate parameters in these fits hence the 
large absolute values of reflection fraction quoted in Table 
[I] The relatively poor chi-squared value of the model quoted 
in Table [l] is due to the higher energy HEXTE data. Re- 
moving the HEXTE data and fitting only PCA data results 
in a x 2 /d.o.f of 16.56/18, but it was considered important 
to include the HEXTE data during fitting to better model 
the electron cut-off and thermal Comptonisation continuum. 
Grouping the HEXTE data in to bins of 4 channels each in 
this fit improves the xVdof to 38.44/28, a perfectly accept- 
able fit given the low number of degrees of freedom. Previous 
pulse-profile modeling ( Poutanen & Gierliriski 2003 1 has de- 
rived a lower limit of 65° for the inclination of this source, 
but since no eclipses or absorption dips are seen, a value of 
60° is assumed throughout this work consistent with pre- 
vious analysis (Gicrlihski, Done, & Barret 2002, Ibragimov] 



& Poutanen 20091. Any difference in inclination only in- 



troduced a systematic shift in the normalisation of spectral 
components, and therefore did not effect the relative phases 
of any pulsations. 
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Figure 1. Top: Unfolded phase-averaged spectrum (model 
multiplied by the ratio of observed data to the re- 
sponse folded spectral model) showing iron line, black- 
body and reflected power-law components using the 
consxphabsx (nthcomp+bbodyrad+rdblurx (pexriv+gaussian)) 
model. Bottom: Data to model ratio plot of this fit. 



3.2 Phase- resolved energy spectra 

3.2.1 Evidence for a varying iron-line EW with phase 

A simple phenomenological fit was performed to 
the 16-bin phase-resolved data consisting of a 
phabs*(pow+bbodyrad) model. The absorption was 
fixed at 0.113 x 



10 cm as in Ibragimov & Poutanen 
(2009) (obtained from the HEADAS tool nh) and the fit 
performed over the region 3.0-25.0 keV whilst ignoring the 
iron-line region of 5.5-10.0 keV. Component normalisations 
were allowed to be free for each of the 16 different spectra. 
Having fitted a model excluding the data in the iron-line 
region, by then including the data in this region one can 
plot the data-to-model residuals as in Figure [2] The plot 
of the residuals demonstrates a broad iron-line feature and 
hints at a changing iron-line equivalent width with phase. 

In the following sections, we shall explore different mod- 
els to explain this apparent pulsation of the iron-line with 
phase. 



3.2.2 Phenomenological Fits 

A fit was performed with the model 
phabs*(cutoffpl+bbodyrad+diskline), using the 
16 phase-resolved spectra over 3.0 to 25.0 keV and the 
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Figure 2. Data to phabs*(pOW+bbodyrad) model residuals 
showing a broad excess in the iron line region. For clarity only 
three of the sixteen phase bins are shown here. The slight soft ex- 
cess is most likely due to the Wien tail of accretion disc emission. 
The zero point from PCU channel 10 has been removed from this 
plot as discussed in the last paragraph of section [2] 



phase-averaged HEXTE spectrum over 15.0 to 150.0 keV. 
The iron line energy was constrained to lie within 6.4 and 
6.9 keV and the inclination fixed at 60° as in Gicrlifiski, 



Done, & Barret (20021 and Ibragimov & Poutanen (2009) 



With the cut-off power-law and blackbody normalisations 
free to vary between phase-bins (untied), and the diskline 
normalisation constrained to its best fitting value across 
all phase bins (tied), the fit gave a x 2 /dof of 756.52/467. 
The best fitting diskline R in was 11.9 R g and the low 
value (1.63) of power-law photon index (compared to the 
phase-averaged fits) points to some degeneracy between the 
power-law and blackbody contributions at soft energies in 
this simple fit. Untying the diskline normalisation between 
phase-bins improved the x 2 /dof to 694.46/451 and the 
spectral components demonstrate the pulsations shown 
in Figure [3] Untying the diskline normalisation simply 
allows the flux from the line to change with phase, rather 
than assuming the flux is constant as a function of phase. 
If the blackbody temperature is untied between phase-bins 
instead of the diskline normalisation, the x 2 /dof improves 
to 691.38/451 and the spectral components this time 
pulsate as shown in Figure [4] Untying the blackbody 
temperature produces an ~ 80° phase-shift in the pulsation 
of the blackbody normalisation compared to Figure [3j 
suggesting that there is some offset between the Doppler 
shift of blackbody photons (that affects the temperature 
measured by an observer at infinity and is dependent on the 
maximum line of sight velocity) and the projected area of 
the hot-spot which dominates the blackbody normalisation. 
It is also interesting to note that the diskline normalisation 
and kT are in phase in Figures [3] and [4] respectively. Again, 
it should be pointed out that the poor x 2 /Aoi is due to 
the higher energy HEXTE data used to constrain the 
parameters in this fit. If a fit is performed excluding the 
HEXTE data (fixing the now unconstrained high energy 
cut off to the phase-averaged values) the xVdof improves 
to 473.58/407. Although the fits obtained with this simple 
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"180 degrees later the full illuminated 
0.2 0.4 0.6 o.s 1 region can be seen. 



Figure 3. Phenomenological fit demonstrating the pulsation in 
the DISKLINE normalisation with phase. The blackbody temper- 
ature was tied between phase bins in this fit. The approximate 
peak-to-peak change in the line normalisation was 42%, in the 
blackbody normalisation 20% and in the power-law normalisa- 
tion 11%. 



Figure 5. Shielding of the reflected emission by the neutron star 
itself could lead to a variation in the visible area of line emission 
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Figure 4. Phenomenological fit demonstrating the pulsation in 
the blackbody temperature when DISKLINE normalisation is tied 
between phase-bins. The approximate peak-to-peak change in the 
blackbody temperature was 5%, in the blackbody normalisation 
28% and in the power-law normalisation 11% 



phenomenological model are formally unacceptable, they 
do provide some physical insight in to the nature of the 
pulsations in this source which are discussed below. The 
detailed reflection models explored in section [3.2.3| produce 
much more acceptable values of reduced \ 2 - 

This preliminary analysis leads to two explanations for 
the spectral pulsation, which are equally likely on statisti- 
cal grounds alone. In the first case, the line normalisation 
is pulsating due to solid angle changes of the line emitting 
region and in the second case the line normalisation itself 
is constant, but the changing spectral shape of the black- 
body component (governed by temperature changes, possi- 
bly due to Doppler effects) introduces a modulation at the 
spin frequency. It is important to point out here the phys- 
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Illuminated disc effectively unchanged with phase. 

Figure 6. Illustration of a geometry that could give rise to a 
constant reflection and iron line emission with phase. The small 
angle between the shock and the rotation axis means that the 
shock docs not preferentially illuminate parts of the disc with 
different line of sight velocities. The illuminated area of the inner 
disc is illustrated using darker shading. 



ical situation that these alternative models are describing. 
In the first case, which invokes a varying line normalisation, 
the visible area of line emission changes with phase. This 
could, for instance, be due to the neutron star self-shielding 
the line emission from the disc as the system rotates, as 
shown in Figure|5j or due to solid angle changes of the emit- 
ting area with phase. In the second scenario, where the line 
emission is constant with phase, the alignment of the spin 
and magnetic poles are such that the illumination of the 
disc by the accretion shock is more or less constant with 
phase as depicted in Figure [6] The amount of observed iron- 
line emission is therefore fairly constant with phase, and the 
modulation can be explained solely in terms of the spectral 
changes in the blackbody component occurring in the same 
energy channels as the iron-line. 
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3.2.3 Fits with reflection models 



Having established these two alternative scenarios using the 
phenomenological fits, the same physically motivated re- 
flection model described in the phase-averaged fits was fit- 
ted to the phase-resolved spectra (and HEXTE spectra). 
In the first fit, the normalisations of the bbodyrad, nth- 
COMP, pexriv and Gaussian components were untied be- 
tween phase-bins, whereas the blackbody temperature and 
seed photon temperature were tied between phase-bins (and 
tied to each other). Throughout these phase-resolved fits, 
all parameters constrained to not vary with phase remained 
very close to the phase-averaged values quoted in Table [l] 
The x 2 /dof for this fit was 455.0/435 and the pulsation of 
the spectral components with phase can be seen in Figure 

m 

It is clear from Figure [7] that, in this model, the gaus- 
sian (iron-line) normalisation is systematically varying with 
phase, whilst the reflection normalisation shows no clear ev- 
idence of doing so. We jointly fitted both the GAUSSIAN and 
pexriv normalisations shown in Figure[7]with a simple con- 
stant offset plus a sinusoid, which forced the fractional am- 
plitude and phase of the normalisation variations to be the 
same. The fit was poor, but was significantly improved by 
untying the phase (Ax 2 = 3.23 for one additional free pa- 
rameter) and amplitudes (Ax 2 = 4.69 for one additional free 
parameter). The likelihood of an observed improvement by 
chance was 0.003 and 3.7 x 10 -6 for untying phase and then 
fractional amplitude respectively. The fractional amplitude 
of variation in the pexriv normalisation was close to zero 
i.e. consistent with a constant pexriv reflection component. 
This inconsistency in the reflection components is the most 
compelling evidence that the "iron line pulsation" is actually 
no more than a change in the underlying blackbody spectral 
shape with phase. The reflection and iron-line normalisation 
should be strongly correlated in amplitude and in phase with 
one another because they both emerge as a consequence of 
hard X-ray irradiation of the disc. In fact, Figure [7] demon- 
strates no discernible pulsation in the pexriv component at 
all. 

Untying the blackbody temperature between phase-bins 
in this fit did not significantly improve the fit. 

A second fit was performed to the phase-resolved spec- 
tra using the same model described in the phase-averaged 
fits, but this time constraining the reflection and line nor- 
malisations to be constant with phase. The y 2 /dof of this 
fit was 492.46/467. Allowing the blackbody temperature to 
also vary with phase (it was still assumed that the blackbody 
temperature and seed photon temperature for Comptonisa- 
tion were identical) improved the x 2 /dof to 464.07/451 with 
an f-test probability of 0.039. The variations with phase are 
shown in Figure [§] Again, as with the earlier phenomeno- 
logical fits, the blackbody normalisation is shifted along in 
phase when the blackbody/seed photon temperature is a 
free parameter. Figure [9] shows 99.9% confidence contours 
for three phase bins (1,6 and 11) of blackbody temperature 
against blackbody normalisation and demonstrates that the 
temperature and normalisation are correlated. However, the 
fact that the 99.9% contours do not overlap significantly 
in temperature for a given normalisation (and vice versa) 
and there is no anti-correlation in the pulsations of these 
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Figure 7. Pulsations in the normalisations of spectral compo- 
nents using the physical reflection model described in the phase- 
averaged fits. Vertical axes are plotted to roughly show the rela- 
tive amplitude of pulsations and the iron-line is modelled using a 
Gaussian component. 
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Figure 8. Pulsations in the BBODYRAD and NTHCOMP normalisa- 
tions and the blackbody temperature with phase with all reflec- 
tion parameters fixed between phase-bins. 



components in figure [8] suggests that the pulsation in these 
components might be real. 

Untying the seed photon temperature for thermal 
Comptonisation and the blackbody temperature had no sig- 
nificant effect on the spectral fits. 

There is no concrete statistical evidence to favour black- 
body temperature variations over varying reflection, but 
there is a good physical reason to do so as outlined above. 
The best fit arises by invoking a changing blackbody temper- 
ature with phase and, although the error bars are large, the 
amplitude and shape of reflection variations do not match 
the line variations seen in Figure [7] In the next section we 
will show that the observed variations can be explained 
quite simply in terms of the Doppler shifts (to produce 
the observed hotspot temperature variation) and solid angle 
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Figure 9. 99.9% confidence contours of blackbody normalisation 
versus kT corresponding to phase bins 1, 6 and 11 in figure|8] Due 
to the coarseness of the grid used to search parameter space (a 
step size of 0.02 keV and 5 for kT and normalisation respectively) 
the contours do not join up. 



Table 2. Table of fractional rms values for different model com- 
ponents 



Fractional RMS 




Fig 3 


Fig4 


Fig7 


Fig8 


Diskline 


0.157 


n/a 


n/a 


n/a 


Bbody norm 


0.072 


0.10 


0.083 


0.046 


Power-law 


0.039 


0.039 


n/a 


n/a 


kT 


n/a 


0.018 


n/a 


0.015 


nthcomp norm 


n/a 


n/a 


0.044 


0.030 


Pexriv norm 


n/a 


n/a 


0.020 


n/a 


Gauss norm 


n/a 


n/a 


0.155 


n/a 



Fractional RMS values are obtained from the best-fitting 
constant offset plus sinusoid model. Quoted value for the Pexriv 
normalisation is an upper limit. 



changes (to produce the normalisation variation) which are 
expected for this system. 

Table [2] shows the fractional rms values for the different 
model components of each model described above. 



4 DISCUSSION 

The rms values shown in Table[2]for the physically motivated 
reflection model invoking blackbody temperature variations 
(Figure |8j, suggest a larger variation in the blackbody com- 
ponent than the Comptonised component. As shown in Fig- 
ure [ToJ this is most likely due to blackbody emission arising 
from a flatter, more horizontally extended 'pancake' type 
region compared to the Comptonised emission from the ac- 
cretion shock which could extend vertically from the neu- 
tron star surface. As the star rotates, the projected area of 
the blackbody region changes more than the accretion shock 




Blackbody hotspot 



Figure 10. The projected area of the blackbody hotspot changes 
more than the projected area of the accretion shock region as 
the neutron star rotates. A snapshot at three different rotational 
phases arc included on this diagram. 



region, leading to a greater amplitude of variation. The frac- 
tional variation from the Comptonised emission might also 
be reduced if the accretion column is relatively optically 
thin (which mitigates the effect of solid angle variation), 
or if there exists a second, constant Comptonising region, 
e.g. an extended corona. The pulsation in the Comptonised 
component is significantly out of phase with the blackbody 



component, as seen in previous work ( Gierlinski, Done, & 



I Barret 2002 Ibragimov & Poutanen 20091. However, the 



phase shift of ~170° between these components in Figure 
[8] (when kT is a free parameter), is much larger than in 
previous work where the shift was ~50° (Gierlinski, Done, 
|fc Barret|[2002l> and ~70° (|Ibragimov fePoutanen| 120091 



with kT fixed. By fixing the blackbody temperature, the 
offset between the blackbody and Comptonised component 
is found to be ~ 50°, consistent with Gierlins ki, Done, fc| 



Barret (20021. It would therefore appear that some of the 
role of changing normalisation is subsumed by the chang- 
ing kT, shifting the variations in blackbody normalisation 
to later phases. It is perhaps surprising how sensitive the 
phase-shift of the blackbody component pulsations are to 
the choice of model used. If the accretion column is opti- 
cally thick and vertically elongated above the neutron star 
surface (as argued above), then the solid angle of the column 
could be maximised when it points away from the observer. 
This would coincide with the minimum projected hotspot 
area and can therefore explain the ~180° phase shift be- 
tween the blackbody and Comptonised components. 

For a 15 km neutron star, the Doppler shift AE/E— 
2nRv/c, produces a large shift at the equator (e.g., around 
13% at a rotation frequency of 401Hz). For an accretion 
shock aligned at 10° to the rotation axis, the upper limit 



obtained by Ibragimov & Poutanen (2009), R is reduced 



from 15 km at the equator to just 2.6 km and the Doppler 
shift reduced to just 2%. If we discount the possibility of 
pulsed iron line variations (as implied by the much weaker 
variations in the reflection continuum), the phase-resolved 
fits suggest that some of the pulsed variation in the spec- 
tral shape can be explained by a varying Doppler shift of 
the blackbody photons underlying the iron-line. If we as- 
sociate the blackbody temperature change with a Doppler 
shift, Figure [8] demonstrates a peak-to-peak change in tem- 
perature of approximately 0.035 keV which in this source, 
assuming an inclination of 60° , translates to a radius of ro- 
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tation of 3.5 km, consistent with the arguments of |Poutanen| 
& Gierlihski (20031. This shift is consistent with a hot-spot 
aligned at approximately 13.5° (20°) to the rotation axis, 
for a 15 km (10 km) radius neutron star. The fact that the 
shock makes a small angle with the rotation axis of the neu- 
tron star implies that the whole inner disc is illuminated 
more or less evenly i.e. there is no preferential illumination 
of different parts of the disc with different line of sight ve- 
locities which would lead to a change in the iron-line shape 
with phase. 

A phase shift measured as ~ 80° can be seen in Fig- 
ure [5] between the pulsation in the blackbody temperature 
(driven by Doppler variations) and the blackbody normal- 
isation (driven by projected area effects). As dicussed in 



Gierlihski, Done, & Barret (20021, this is consistent with a 



geometry where Doppler variations are strongest when the 
hot-spot is at the limb of the star, and projected area effects 
are strongest ~ 90° later when the solid angle of the hot- 
spot is maximised. Clearly, this picture is only approximate, 
since we do not take account of effects such as light-bending 
(e.g. Ibragimov fc Poutanen|2009 |. 

Our fits to these data are consistent with the idea that 
the iron-line shape and normalisation remain constant with 
phase, suggesting a geometry where the illumination of the 
disc by the accretion shock is fairly constant. By assuming 
that the only variable components are the normalisations 
of the blackbody and Comptonised components (|Gierlihski,| 



Done, & Barret 2002 Ibragimov & Poutanen 2009), one 



can detect an apparent pulsation in the residual iron-line 
component. However, any apparent change in the iron-line 
normalisation with phase can be explained by taking into 
account the Doppler boosting of the underlying blackbody 
continuum. The constancy of the iron line is not really sur- 
prising, because for the small offset of the accretion column 
from the rotation axis that is envisaged here, the disc will 
not see much variation in solid angle of the accretion-shock 
region, and variations in Doppler beaming of the accretion- 
shock emission toward the disc will be small. The presence of 
a constant power-law continuum component from a corona 
would also serve to reduce the variability of the iron line 
emission. 



5 CONCLUSIONS 

Phase-resolved spectroscopy of the 2002 outburst of the 
AMXP SAX J 1808.4-3658 as observed by RXTE reveals no 
convincing evidence for pulsations in the iron-line or changes 
of any line properties with phase. By allowing the blackbody 
spectral shape to change (i.e. freeing the blackbody tem- 
perature during fitting) one can allow for variable Doppler 
boosting of the underlying blackbody continuum, which in- 
troduces a larger phase-shift between the Comptonised com- 
ponent normalisation and the blackbody normalisation than 
has been seen in previous work. The phase-shift between the 
blackbody normalisation and the blackbody temperature is 
approximately consistent with the expected 90° offset be- 
tween the maximum projected area and maximum Doppler 
shifts. The fraction of the disc illuminated by the Compton- 
ising region appears to be roughly constant with rotational 
phase. 
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